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Abstract. We present a general Ginzburg-Landau theory of electrostatic interac- 
tions and electric field effects for the order parameter, the polarization, and the 
charge density. Electric field effects are then investigated in near-critical fluids and 
liquid crystals. Some new predictions are given on the liquid-liquid phase transition 
in polar binary mixtures with a small fraction of ions and the deformation of the 
liquid crystal order around a charged particle. 

1 Introduction 

There can be a variety of electric field effects at phase transitions. For near-critical 
fluids a number of such effects have been investigated [1-19]. Examples are weak 
critical singularity in the macroscopic dielectric constant and the refractive index, 
nonlinear dielectric effects, and electric birefringence. In liquid crystals electric 
field is coupled to the nematic orientation [20], but similar effects exist near the 
isotropic- nematic transition [21]. For near-critical fluids without charges we here- 
after summarize previous results in this section. For other systems such experments 
are not abundant, so we will try to predict some effects in the following sections. 

In near-critical fluids the order parameter if> has been taken to be n/n c — 1 for pure 
(one-component) fluids and to c — c c for binary fluid mixtures, where n is the density 
and c is the molar concentration or volume fraction of one component. The quantities 
with the subscript c will represent the critical values. To be precise, however, there 
is a mapping relationship between fluids and Ising systems in describing the critical 
phenomena [22]. The number density in fluids may be expressed as 

n{r) =n + a 1 ^{r)+(3 1 ^{r) 2 , (1.1) 

in terms of the spin ip(=the true order parameter) in the corresponding Ising system. 
Although not well justified, we expand e in powers of ip in a local form as 

e(r) = e + eiV>(r) + ^e 2 ^(r) 2 + ■■■, (1.2) 



"NATO series Nonlinear Dielectric Phenomena in Complex Liquids (Kluwer, 2003) 



1 



where S\ and 62 are constants independent of the reduced temperature r = T/T c — 1. 
The quantity ip 2 corresponds to the energy density in Ising systems and the average 
(if) 2 ) over the thermal fluctuations contains a term proportional to |T| 1_a on the 
critical isochore above T c and on the coexistence curve below T c [22] . Therefore, the 
weak singularity (oc |r| 1_a ) appears in the so-called coexistence curve diameter for 
Pi 7^ below T c and in the thermal average (e) for e 2 7^ above T c . 

Regarding the expansion (1.2) we make two remarks, (i) For nonpolar pure 
fluids the overall density-dependence of e may well be approximated by the Clausius- 
Mossotti relation (e — + 2) = a m n/3, where a m is the molecular polarizability. 
Using this relation and assuming p 1 — in (1.1) we have Si = (e — l)(e + 2)/3 
and £ 2 = 2ef/(2 + e ) > 0. However, dielectric formulas like Clausius-Mossotti 
are valid only in the long wavelength limit where the electric field fluctuations are 
neglected (in text books) or averaged out [3, 4]. Thus estimation of 62 from Clausius- 
Mossotti is not justified when we use (1.2) to derive dielectric anomaly due to the 
small-scale critical fluctuations. Notice that (1.1) itself gives rise to a contribution 
(= 2/3 1 £ 1 /aii) to e 2 from the linear density dependence e = e c + Si{n/n c — 1). (ii) 
For binary mixtures A+B, e\ should be of order Ea — £b, where e& and e-q are the 
dielectric constants of the two components, so E\ can be 10-100 in polar binary 
mixtures (where at least one component is polar). Electric field effects are much 
stronger in polar binary mixtures than in nonpolar fluids. 

Once we assume (1.2), it is straightforward to solve the Maxwell equations within 
the fluid if the deviation 5e = e — (e) is treated as a small perturbation. Hereafter 
(• • •) represents the average over the fluctuations of tp. The macroscopic static 
dielectric tensor e e n is obtained from the relation (D) =£ c ff ~(E) between the 
average electric induction (D) and the average electric field (E). To leading order, 
the fluctuation contribution Ae eS in the long wavelength limit can be expressed in 
terms of the structure factor I(q) = (\ipq\ 2 ) [3, 23]. The contribution from the long 
wavelength fluctuations with q < A is written as 



where jq = (2tt)~ 3 f dq, 1= {5^} is the unit tensor, and q = q~ 1 q. The upper 
cutoff wave number A in (1.3) is smaller than the inverse particle size. Note that 
I(q) is uniaxial at small q in the direction of electric field, but it is nearly isotropic 

at large q. Thus the contribution Ae eS from the short wavelength fluctuations with 
q > A is proportional to the unit tensor and may be included into eo in (1.2). 

We may also calculate an electromagnetic wave within near-critical fluids [3, 23]. 
Let k be its wave vector. The average electric field (E) is nearly perpendicular to 
k, oscillates with frequency u, and obeys 




(1.3) 



oo 2 (e(w)) + Ae(fc,w) ■ (E) = c 2 k 2 (E), 



(1.4) 
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where e(uS) = e (uj) + ei(uj)ip + ••• is the dielectric constant at frequency u(= 
ckj Eq(uj) 1 I 2 ) dependent on ip as in (1.2). The fluctuation contribution to the dielec- 
tric tensor for the electromagnetic waves is written as [3, 23, 24] 

A7(k,.) = ^f J^pLjk^,-^,, (1.5) 
e (u) Jq q 2 — k 2 — tO 

where q ± = q — k~ 2 (q ■ k)k and I ±— I — k~ 2 kk are perpendicular to k, and +i0 
represents a small imaginary part arising from the causality law. The imaginary 
part gives rise to damping and can be calculated using the formula 1/ (x — iO) = 
vp(l/x) + i7i5(x) with vp representing taking the Cauchy principal value. Then 
Ae(k,u) has two eigenvalues A a (co) (a = 1,2) in the plane perpendicular to k. 
The dispersion relations for the two polarizations are written as 

k = c- 1 u[(e(uj)} + A ± (uj)} 1 / 2 . (1.6) 

In particular, if we apply an electric field along the z axis and send a laser beam 
along the x axis, we have Ae yz = Ae zy = 0. The principal polarization is then either 
along the z axis (a = 1) or along the y axis (a = 2). The turbidity T a is written as 

T a = — ^-rImA a (w) = -^T^iH 2 / dQI(k - kq)(l - q 2 a ), (1.7) 

where q\ = q z , <j2 = q y , and dQ represents integration over the direction q. The 
Ornstein-Zernike intensity gives the famous turbidity expression. Generally, when 
I(q) depends on q, the difference Ai — A 2 becomes nonvanishing. Its real and 
imaginary parts can be measured as (form) birefringence and dichroism, respectively. 
In near-critical fluids under shear flow, both these effects exhibit strong critical 
divergence [25]. In polymer science, birefringence in shear flow arising from the above 
origin is called form birefringence [24], whereas alignment of optically anisotropic 
particles gives rise to anisotropy in (e^-) leading to intrinsic birefringence. Form 
dichroism is maximum when the scattering objects have sizes of the order of the 
laser wavelength and, for shear flow, it has been used to detect anisotropy of the 
particle pair correlations in colloidal suspensions [26] and polymer solutions [27]. In 
Ref. [23] an experimental setup to measure dichroism was illustrated. 

In Section 3 we will examine the effects mentioned above using (1.3) and (1.5). 

2 Ginzburg-Landau free energy in electric field 

We will construct a Ginzburg-Landau free energy including electric field under each 
given boundary condition, since there seems to be no clear-cut argument on its form 
in the literature. In our scheme the gross variables are the order parameter ip, the 
polarization p, and the charge density p. They are coarse-grained variables with 
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Figure 1: (a) System of a capacitor and a dielectric material with inhomogeneous dielectric 
constant e(r) at fixed capacitor charge Q. The potential difference $ is a fluctuating 
quantity dependent on e(r). (b) Two capacitors connected in parallel with charges Q and 
Qq. The smaller one contains an inhomogeneous dielectric material, and the larger one 
a homogeneous dielectric material. In the limit Q/Qo — > 0, the potential difference 
becomes fixed, while Q is a fluctuating quantity. 

their microscopic space variations being smoothed out. Since the electromagnetic 
field is determined for instantaneous values of the gross variables, we will suppress 
their time dependence. 

2.1 Dielectrics under given capacitor charge 

The first typical experimental geometry is shown in Fig. la, where we insert our 
system between two parallel metallic plates with area S and separation distance L. 
We assume S 1 ^ 2 ^> L and neglect the effects of edge fields. Generalization to other 
geometries is straightforward. The z axis is taken perpendicularly to the plates. Let 
the average surface charge density of the upper plate be a ex and that of the lower 
plate be — cr ex . The total charge on the upper plate is 

Q = Sa ex . (2.1) 

The Ginzburg-Landau free energy functional F consists of a chemical part F c ^{ip, p, p] 
and an electrostatic part F st {p, p,Q} as 

F = F ch {^,p,p} + F st {p,p,Q}, (2.2) 

where if) represents a set of variables including the order parameter. In ferroelectric 
systems p is the order parameter. The equilibrium distribution of the gross vari- 
ables is given by const. exp(—F/k-BT) at each fixed Q. We determine F st as follows. 
If infinitesimal deviations 5p, dp, and SQ are superimposed on p, p, and Q, the 
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incremental change of F st should be given by the work done by the electric field, 

5F st = F st {p + 5p,p + 5p,Q + 5Q} - F st {p,p,Q} 

= J dr[-E ■ 5p + (f)5p} + <S>5Q, (2.3) 

where the space integral is within the system between the plates, E = — V(/> is the 
electric field vector, and <3> is the potential difference between the two plates. The 
electric potential may be set equal to at the lower plate and $ at the upper 
plate. The electric induction D = E + AnP satisfies 

V D = -V 2 + 4ttV • p = 4vrp, (2.4) 

in the bulk region. The potential satisfies 

-V 2 = 47rp eff , (2.5) 

where 

p cff (r) = p{r) - V • p(r) (2.6) 
is the effective charge density. The boundary conditions at z = and L are 

E x = E y = 0, D z = -47r<7 ex . (2.7) 

With these relations of electrostatics we can integrate (2.3) formally as 

= ±fdr#. (2.8) 



E 



St 



In fact (2.8) leads to the second line of (2.3) if use is made of S(E 2 ) = -2V0 • 5D - 
8%E ■ 5p = -2V((f)SD) + 8n((p5p -E-Sp). 

To explicitly express F st in terms of the gross variables, we first assume that all 
the quantities depend only on z and p is along the z axis, for simplicity [28] . If we 
define w(z) = Jq dz'p(z'), we obtain 

4>{z) = 4tt f dz'\p z (z') + a cx - w(z% (2.9) 
Jo 

From the overall charge neutrality condition we require w(L) = 0, so (2.7) is satisfied 
and 

F st = ^ / dzE z {zf = 2rrS [ dz[ Pz {z) + a cx - w{z)} 2 . (2.10) 

87T JO JO 

For general inhomogeneous p and p, we define the lateral averages, 

P(z) = gj dr ± p(r), p(z) = -j dr ± p(r), (2.11) 
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where r± = (x, y) is the lateral position vector. We may assume p x = p y = from 
the geometrical symmetry in the limit S 1 ^ 2 / L — > oo. The effective charge density 
p c s is divided into the lateral average p — dp z /dz and the inhomogeneous part, 

Pinh(r-) = W) - p{z)] - V • \p(r) - p(z)}. (2.12) 

The electric potential is expressed as 

<t>{r) = An £ dz'[p z (z') + a cx -w(z')}+An J dr'G(r,r')p inh (r'). (2.13) 

The first term is of the same form as (2.9) and 

w(z) = f dz'p(z'), (2.14) 



which vanishes at z = and L. The Green function G(r, r') = G(r', r) satisfies 

V 2 G(r, r') = -5(r - r') (2.15) 
and vanishes as r and r' approach the surfaces of the conductors as 

lim G(r,r')= lim G(r, r') = 0. (2.16) 

Then we have E x = E y = at z = and L. The potential difference $ is written as 

$ = 4tt [ dz[p z (z) + a cx - w(z)\ 
Jo 

= 4TvL(p z + a cx + zp) s , (2.17) 

where (• • -) s = / dr(- ■ •)/ SL represents the space average and use has been made of 
(p) s = 0. This relation may also be written as 

$/4ttL - (p z ) s = Q/S - ((L/2 - z)p) s , (2.18) 

which relates the mean electric field — $/L and the capacitor charge Q. The second 
term on the right hand side becomes important when the charge is accumulated 
near the capacitor plates in electric field. 

Because we are taking the limit S 1 ^ 2 / L — > oo, the translational invariance 
G(r,r') = G(rj_ — r' ± ,z,z') holds in the xy plane. The Fourier transformation 
Gk(z, z') in the xy plane satisfies (—V 2 + k 2 )Gk(z, z') = 5(z — z') and becomes [8] 



G k (z,z') = G k (z>,z) = ^-e~ k ^~ 



x 



2k 2ksinh(kL) 

smh(kz)e- k{L - z>) + smh(kL - kz)e~ kz ' 



(2.19) 
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In particular, in the long wavelength limit k — > 0, we have 



G (z,z') = lim G k (z,z') = -[z + z' - \z — z'\\ — zz'/L. 



(2.20) 



From (2.17) and (2.20), as it should be the case, the potential is also expressed as 

(f)(r) = <$>z/L + 4tt J dr'G(r, r')p eS {r'). (2.21) 
The electrostatic energy is written in terms of p in h as 



F st = 2nS I dz 



Jo 



p z (z) + a ex -w(z) 



+ 



2tt jdrj dr'p mXl (r)G(r,r')p mXl (r'). (2.22) 



With the aid of (2.17) and (2.20), F st in terms of p cS reads 

S 



F st = — $ 2 + 2tt Jdrj dr'p cS (r)G(r, r')p cS (r'). (2.23) 



2.2 Dielectrics under given capacitor potential 

In the previous case, the capacitor charge Q is a control parameter and the potential 
difference $ is a fluctuating quantity. We may also control $ by using (i) a battery at 
a fixed potential difference connected to the capacitor or (ii) another large capacitor 
connected in parallel to the capacitor containing a dielectric material as in Fig. lb. 
We examine the well-defined case (ii). The area So and the charge Qq of the large 
capacitor are much larger than S and Q, respectively, so the large capacitor acts 
as a charge reservoir. We are supposing an experiment in which the total charge 
Qtot = Qo + Q is fixed and the potential difference is commonly given by $ = 
Qo/Cq, where Cq is the capacitance of the large capacitor. Obviously, in the limit 
Q/Qo ~ S/So — > 0, the deviation of $ from the upper bound $tot = Qtot/ Co 
becomes negligible. Because the electrostatic energy of the large capacitor is given 
by E = Ql/2C = (Q t 2 ot /2C )(1 - Q/Q tot ) 2 , we obtain 

E = (Q 2 tot /2Co) - $Q, (2.24) 

where the first term is constant and a term of order (Q/Qtot) 2 is neglected. There- 
fore, for the total system including the two capacitors, the potential difference $ 
is a control parameter and the appropriate free energy functional is given by the 
Legendre transformation [1], 

G = F-$Q. (2.25) 

If $ is fixed, Q is a fluctuating quantity determined by (2.18). The equilibrium 
distribution of the gross variables is given by const. exp(— G/A; B T) at each fixed $. 
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2.3 Dielectrics in vacuum and dipolar ferromagnets 

There is another typical case, in which a sample of dielectrics is placed in vacuum. 
There is no polarization and charge immediately outside the sample, but there can 
be an applied electric field created by capacitor plates far from the sample. As far 
as I am aware, simulations of phase transitions in dipolar particle systems have been 
performed under this boundary condition [29], where the particles are assumed to 
interact via the dipolar potential oc Sij/r 3 — 3xiXj/r 5 even close to the surface. In 
ferromagnetic systems there is no magnetic charge and this boundary condition is 
always assumed [30]. 

We assume no free charge (p = 0) and apply an electric field E ex along the z 
axis. The electrostatic potential is given by 

<j>(r) = -E cx z + J dr f G (r,r')p cS (r r ) + J daG (r,r a )a cfi (r a ) 

= -E cx z + J dr'p(r') ■ VG (r - r'). (2.26) 

In the first line p eff = —V • p is the effective charge as given in (2.6), a cS = p ■ n 
is the effective surface charge with n being the outward normal unit vector at the 
surface, and r a represents the surface point. In the second line the space integral 
is within the sample and V' = d/dr'. Since there is no capacitor charge in contact 
with the system, Go is of the usual Coulombic form, 

Go(r-r') = l/\r-r'\. (2.27) 

The total free energy consists of the chemical part and electrostatic part as F = 
F c h{^P,p} + Ft{p}- The latter is the usual dipolar interaction, 

F Bt = -J drE cxPz + \j dr j dr' ^Pi^pM')^ N jG (r - r% (2.28) 

ij 

which satisfies SF st /5p = —E in (2.3). Here — VjVjG^r) = <W rS — SxiXj/r 5 , 
as ought to be the case. In the ferromagnetic case, p and E cx correspond to the 
magnetization and the applied magnetic field, respectively [1]. 

The above electrostatic (magnetic) energy (2.28) depends on the sample shape 
due to the depolarization (demagnetization) effect, so the phase transition to a 
ferroelectric (ferromagnetic) phase becomes shape-dependent [29] . As is well known, 
if the dielectric constant e is uniform in the sample, the field inside is given by E ex /e 
for a thin plate and 3E CX / (2 + e) for a sphere along the z axis. 

3 Near-critical fluids without ions 

As an application of the general scheme in Section 2, we consider a near-critical fluid 
in the geometry of Fig. la. We will use the results already presented in Section 1. 
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3.1 General relations in the charge-free case 

In the absence of free charges (p = 0), F ch in (2.2) is written as 

F ch {iP,p} = F {iP} + J dr^-p 2 , (3.1) 

where X represents the electric susceptibility. The first term is the usual free energy 
functional for the order parameter ip. Close to the critical point it is of the form, 



FoW = k B T c J dr 



—rip + —ip +— \Vip\ 



(3.2) 



where a , u and K are positive constants, and r is the reduced temperature. The 
term oc / drip is not written (because ip is a conserved variable). Use of (2.3) gives 



6F = 5F + J dr [(rt -E)-5p- ^p 2 S X 



+ $6Q. (3.3) 



2x 2 

Thus F is minimized for 

p = X E, D = eE. (3.4) 
The static dielectric constant is 

e = l + 47vx (3.5) 

and depends locally on the order parameter ip as in (1.2). In usual fluids p rapidly 
relaxes to xE on a microscopic time scale, so (3.4) may be assumed even in nonequi- 
librium. Then p is determined by ip and Q. The charge- free condition within the 
fluid V • D = is written as 

V • (eV0) = 0. (3.6) 
From the general expressions (2.8) and (3.1) with the aid of (3.4) we obtain 

F = F {iP} + j-jdrE-D = F {ip} + i$Q. (3.7) 

The expression in terms of Q also follows directly from integration of 5F = &5Q at 
fixed ip, where the factor 1/2 arises because the ratio <3>/Q is a functional of ip only 
and is independent of Q from (3.4). On the other hand, in the fixed-potential case 
in the geometry of Fig. lb, the free energy functional G in (2.25) is written as 

G = F {iP} I drE-D = F {iP} - Uq. (3.8) 

on J Z 

The electrostatic parts AF = F — F and AG = G — F are functionals of e(r) at 
fixed Q and $, respectively. The functional derivative of AF at fixed Q and that of 
AG at fixed $ with respect to e are both of the same form, 

», = (^). = -^' ^ 
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because / drE -SD=-J drVcp ■ 5D = at fixed Q and / drD ■ 5E = at fixed $. 
These relations directly follow from (3.3) under p = \E. Even if a sample is placed 
in vacuum, (5AF/5e)e assumes the same form. To derive (3.9) more evidently, 
we may assume that e = e(z) is a function of z only. Then we explicitly obtain 
AF = 2nQ 2 dze(z)- 1 and AG = -(Bn^S^ 2 / J L dze(z)- 1 to confirm (3.9). 

3.2 Landau expansion and critical behavior 

We assume that a constant electric field E = (0, 0, — E ) is applied in the negative 
z direction. Then e 2 E 2 /4iTk B T c (= A t c in (3.19)) and (e 2 1 /s )E 2 /A7rk B T c (^ g e in 
(3.16)) are two relevant parameters representing the influence of electric field. We 
assume that they are independent of the upper cut-off A of the coarse-graining, while 
the coefficients in Fq in (3.2) depend on A [22]. This causes some delicate issues. 
Theory should be made such that the observable quantities are independent of A. 

We expand the electrostatic free energy with respect to the order parameter 
ip. In the following we assume the fixed charge condition, but the same expressions 
follow in the fixed potential condition owing to (3.9). Let the electric field be written 
as E = Eq — where 5<p is the deviation of the electric potential induced by 

5e = e — (e). Then (3.9) is expanded as 

To first order in 5e the charge-free condition (3.6) becomes (e)V 2 <50 = E ■ VSe. 
Using the Green function G(r,r') in (2.15) and (e) = e , we obtain 

S(j>(r) = -— *— / dr'G(r, r')E ■ V'fe(r'). (3.11) 

The electrostatic free energy AF = F — F is composed of two parts up to order 
0(5e 2 ). The first part is 

F eo = -^J drE 2 e = ~J drE%(e + + ^ 2 ). (3.12) 

The linear term (oc / drip) is irrelevant for fluids and the bilinear term (oc / drip 2 ) 
gives rise to a shift of the critical temperature as will be shown in (3.19) and (3.22). 
The second part is a dipolar interaction arising from the second term in (3.10), 

F dip = ^A S E 2 J dr J dr'[V z iP(r)]G(r,r')[V z iP(r% (3.13) 
where V z = d/dz. The coefficient A s is defined by 

A s = ej/e = e\/e Q . (3.14) 
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which is small for pure fluids (= (e — l) 2 (^o + 2) 2 /9e from Clausius-Mossotti) , 
but can be 10 — 100 for polar binary mixtures. The F&p is positive-definite for the 
fluctuations varying along Eq but vanishes for those varying perpendicularly to Eq. 
Thus F dip produces no shift of T c in the mean-field theory, but its suppresses the 
fluctuations leading to a fluctuation contribution to the shift as in (3.23) below. 

Far from the capacitor plates we may set G(r,r') = 1/An\r — r'\. Picking up 
only the Fourier components ipq = J drip(r)e l( l' r with qL 3> 1, we obtain 

i^dip = ^k B T c g c J^q 2 z \ij q \ 2 , (3.15) 

where q z = q z /q. The strength of the interaction is given by 

g e = (4ttA;bT c )- 1 A s Eo 2 . (3.16) 

A similar dipolar interaction is well-known for uniaxial ferromagnets [30, 33]. At 
the critical density (or composition) in one-phase states, the structure factor I{q) 
in the presence of Fdi P becomes uniaxial as 

I(q) = [a (r - r c ) + g c q 2 + K q 2 }-\ (3.17) 

where r c is a shift induced by electric field. For g c q 2 > aoT c , the intensity decreases 
with increasing electric field even if r c > 0. 

(i) Weak field regime: If g e is smaller than a^r = x' 1 , the electric field is weak 
and the structure factor may be expanded as 

I(q) = Io(q) + Io(q) 2 (a r c - g c q 2 z ) + • • • , (3.18) 

where I (q) — x/(l + g 2 £ 2 ) is the intensity at zero electric field with x — ^(0)- The 
apparent shift arises from the bilinear term (oc ip 2 ) in F e0 in (3.12): 

a r c = e 2 E 2 /4irk B T c = e 2 g e /A s . (3.19) 

Note that a$T c and hence (3.18) are independent of A. The usual critical behavior 
at long wavelengths g£ < 1 follows for A£ ~ 1, where £ = £+ot v is the correlation 
length. The asymptotic scaling relations are x — Tr -7 (or ao = r _1 r 7_1 ) and 
K = x~ l i 2 = r-^r - ^, where T is a constant, and 7 = 1.24, v = 0.625, and 
rj = 0.03 — 0.05 are the usual critical exponents. 

(ii) Strong field regime: The electric field is strong for r — t c < r c . The crossover 
reduced temperature r e and the characteristic wave number k e are written as 

r e = (Tg c f\ k c = C + l(Tg c yl\ (3.20) 

As a rough estimate, we set Y ~ £j. , £ +0 ~ 2.5 x 10 _8 cm, and T c ~ 300K to 
obtain r e ~ lO -6 ^ 8 ^ 6 with Eq in units of 10 4 V/cm. There is no experimental 
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report in this regime. According to the renormalization group theory in the presence 
of the uniaxial dipolar interaction (3.15) (which treats Fq + i 7 di P )[33], there is no 
renormalization of the coefficient g e and the dipolar interaction strongly suppresses 
the fluctuations with wave numbers smaller than k e in the temperature range r— r c < 
r e resulting in mean field critical behavior for d > d c = 3. For these fluctuations we 
may set A ~ k c to obtain the renormalized coefficients, 

a = r-Vr 1 , i^o^r 1 ^, u = g*K 2 K/K d , (3.21) 

where K d = {2n)~ d 2n d / 2 /T(d/2), e = 4 - d, and g* = e/9 + 0(e 2 ) is the fixed point 
value of the coupling constant g = K d u /KQ A e [22]. The shift r c in (3.17) in strong 
field consists of two terms as r c = r c i + r C 2. The origin of r c i is the same as that of 
r c in (3.19). Here using (3.20) and (3.21) we obtain 

r c i = e 2 El/Ana Q k B T = (e 2 /A s )r e . (3.22) 

The t C 2 arises from the quartic term (oc Mo) in -^o in (3-2): 

3 r 

a Jq 

where C d = K d l £ 2 J£ 2 (£ 4 + £ 2 Z ) (= ?r/4 at d = 3). The r c2 is a positive fluctuation 
contribution meaningful only in strong field. (A similar shift was calculated for 
near-critical fluids in shear flow [22].) Therefore both r cl and r c2 are of order r e . 
We also consider the correlation function g(r) = Jql(q)e 1( l' r in strong field. From 
(3.17) and (3.21) the scaling form g (r) = G(x/£±, y/£±, z/£\\)/K £\\ holds with 

£i = K^/ir ~ r c )] 1/2 , C|| = Ul = K'rJir - r c ) > £±. (3.24) 

Here G(0, 0, Z) ~ Z~ 3 for Z ^> 1. The large-scale critical fluctuations (>> A;" 1 ) are 
elongated along the z axis, as expected in experimental papers [13]. 
We now comment on previous work. 

(i) If do is treated as a constant, T c r cl is the mean field shift in Landau-Lifshitz's 
book for pure fluids [1]. The same shift was later proposed for binary mixtures [7]. 

(ii) Debye and Kleboth [2] derived a reverse shift (AT C = — T c r c i) neglecting the 
inhomogeneity of E and setting E = E in (3.7). The dipolar interaction was 
nonexistent, leading to the normalized turbidity change T(E )/T(0) — 1 = x a o r c 
for !;( < 1, where T(E ) is the turbidity in electric field (see (1.7)). They found 
turbidity decreases in nitrobenzene-|- 2,2,4-trimethylpentane to obtain r c = —0.5 x 
10~ 4 at Eq = 45kV/cm in agreement with their theory, where Hildebrand's theory 
gave do = (4.3A)^ 3 and data of e(<p) as a function of the volume fraction <p of 
nitrobenzene gave d 2 e/d(j) 2 = 28.7. (Notice that their |r c | is of order r e in (3.20) if 
A s ~ 10.) Subsequent light scattering experiments detected similar suppression in a 



1 



1 



OnQ 2 + KnQ 2 KnO 2 



3g*C d r c , 



(3.23) 
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near-critical binary mixture [31] and in a polymer solution [32]. In our theory (3.18) 
and (3.19) hold in weak field, so for fc( < 1 we have 

T(E )/T(0)-l = x(a r c -g c /5) or X (a r c - 2g e /5) (3.25) 

for the polarization along the z or y axis. Thus we predict that the turbidity 
decreases for g c /5 > a r c or for A s > 5e 2 for any polarization. Here we expect that 
A s is considerably larger than \e 2 \ in polar mixtures (see the end of Subsection 3.3 
also). Complex effects also arise from a small amount of ions which are present in 
most binary mixtures (see Section 4). 



3.3 Macroscopic dielectric constant 

@ 

We calculate the macroscopic dielectric constant e e fr = ^ir(Q) / SE = {D z )/E 
in the absence of charges. From the charge-free condition V ■ (D) = W Z (D Z ) = 0, 
(D z ) is homogeneous in the fluid as 

(D z ) = (e)E - (5eV z 5<p) = const. (3.26) 

Far from the boundary surfaces, by setting M = (ip) , we obtain 

I 

'q 



e eS = e + ExM + \e 2 M 2 + f 
2 Jn 



if 2 ~ 



I(q), (3.27) 



where the last term is the zz component of the tensor in (1.3). 

(i) Linear response: The linear dielectric constant £ e ff(0) is defined in the limit 
E — * 0, where q 2 z in (3.27) may be replaced by 1/3. Because a ^ 2 corresponds to the 
energy density in the corresponding Ising model (see the sentences below (1.2)), the 
renormalization group theory gives (aoip 2 ) = const. — At 1 ^ / (1 — a) + • ■ ■ with a = 
0.11 (the specific-heat critical exponent) and A = Rf/(cx^ ) at the critical density 
or composition above T c [22]. In terms of this A, the specific heat Cy for pure fluids 
and Cp C for binary mixtures per unit volume grow as —kBd{aoip 2 )/dr = k-QAT~ a . 
The -R^ = aA^ (= -fQ/4) is a universal number (the two-scale-factor universality). 
It follows that £ e ff(0) should exhibit weakly singular behavior [4-6,8], 

£eff(0) = € c + D 1 r l - a + D 2 r + ---. (3.28) 

In our theory D\ = (A s — 3e2/2)A/3a depends on the arbitrary cut-off A since 
a oc A^ 7-1 ^ for A ^> indicating inadequacy of our theory at short wavelengths. 
To calculate D 1 correctly we need to interpolate the renormalization group theory to 
a microscopic theory [34]. On the other hand, the refractive index n(u) = e(oo) l/2 at 
optical frequency has also been predicted to be of the form of (3.28) [4, 6]. However, 
despite a number of experiments, unambiguous detection of the weak singularity in 
these quantities has been difficult for both pure fluids and binary mixtures [10]. 
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Figure 2: Nonlinear dielectric constant Ae e g(r, Eq) (divided by a constant Bq) vs reduced 
temperature r' = r — t c > in applied electric field Eq for r e = 10~ 5 , 10~ 6 , and 10~ 7 , 
where r e oc Eq® from (3.20). 



(i) Nonlinear response: The nonlinear dielectric constant Ae e s = e e s(E ) — £ e ff(0) 
has been observed to become positive and grow near the critical point in polar binary 
mixtures [12-15]. Because small-scale fluctuations are insensitive to electric field, 
such critical anomaly arises from nonlinear effects at long wavelengths (q <£ _1 for 
weak field and q <k e for strong field). Its calculation is therefore much easier than 
that of £ e ff(0). Use of the structure factor (3.17) yields 

Ae cS -± r (2A B /5-e 2 + el/A s )K- 2 g e £ (r > r c ) (3.29) 

' (A s -e 2 )K Q - 1 k e (0<r'<r e ), (3.30) 



167T 

In the first line we have K oc £ v and Ae e fj oc E 2 ^ l ~ 2r] , while in the second line 
K oc k-* 1 and Ae cS oc E^ n)/<? ~ v) . Here we set r' = r - r c > 0. Then r' = r for 
weak field. If Ae eff is positive for any t', we need to require A s > e 2 . Previously we 
calculated (3.29) without the terms involving e 2 [8]. For general r'/r c we find the 
scaling, 

Ae cS = B r^F non (k e 0, (3.31) 

where B is a constant, -F non (x) ~ x 1_2r? for i < 1, and -F non (a;) — > 1 for x 1. 
Here a; = fc e £ = (t c /t') u . If we neglect the contributions involving e 2 and set r/ = 0, 
F n0 n{x) is calculated in a universal form, 

F non (x) = (I + JL) v^T^ - ^ In(x + v/TT^) - ^. (3.32) 

In Fig. 2 we plot A£ e fj/i?o = 7"g -F nc ,n(&e£) vs t' by neglecting r/. For weak field 
r c <C r, we predict Ae fj oc T~ Zaou E 2 at the critical density or composition, where 
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Znon = v(l ~ 2 v) = 7 - 2 /5 = 0.58 using 7 + 2/3 = 2 - a = Zv [4, 5, 8]. For strong 
field r c >t', Ae cff should saturate into (3.30). 

In the experiments, around 0.4 and the coefficient Ae c s/[r Znon E 2 ] was 

proportional to (e a — Eb) 4 /e 2 for various binary mixtures [13, 14]. The latter aspect 
is in accord with (3.29) if S\ ~ Ea — Eb and A s — 5(e 2 — eI/A s )/2 ~ A s . We also 
note that a weak tendency of saturation of Ae c g was detected at r = 3 x 10~ 4 
with increasing E above 60 kV/cm, but it was attributed to a negative critical 
temperature shift [15]. These results and Debye-Kleboth's data suggest that \e2\ is 
considerably smaller than A s at least in polar mixtures. 

3.4 Critical electrostriction 

We consider equilibrium of pure fluids in which (E) = E (r) varies slowly in space. 
In equilibrium the chemical potential defined by the following is homogeneous [1], 

/i(r, M, E Q ) = (SF/Sn) = /i(t, M, 0) - ^-e^ 2 = const. (3.33) 

Here we set ip = (n — n c )/n c . When E varies slowly compared with the correlation 
length, an inhomogeneous average density variation is induced as 

(Sn)(r) = n c K T £ 1 E (r) 2 /87r, (3.34) 

where Kt = (dn/d^T/n 2 ~ r~ 7 /k^,T c n c is the isothermal compressibility. Experi- 
mentally, the above relation was confirmed optically for SF 6 around a wire conductor 
[35], and was used to determine /i(T, M) for 3 He in a cell within which a parallel- 
plate capacitor was immersed [36]. 

This problem should be of great importance on smaller scales particularly in 
near-critical polar binary mixtures. For example, let us consider a spherical particle 
with charge Ze and radius R placed at the origin of the reference frame. The fluid 
is in a one-phase state and is at the critical density or composition far from the 
particle, so r > 0. In the Ginzburg-Landau scheme (3.2) and (3.9) yield 

(a r + u^ 2 - K Q V 2 )i> = h(r)^ 3 = A Q i +Q r~\ (3.35) 

where A = e 1 Z 2 e 2 /87r£g/cB^+o- The h(r) = A (£ +0 /r) 4 is a dimensionless space- 
dependent ordering field. We may assume E\ > and A > 0. For polar binary 
mixtures and/or for colloidal particles with Z ^> 1, Aq can well exceed 1 (see Section 
4). In the space region, where £ < ((91n/i(r)/(9r) _1 ~ r, the usual scaling relations 
hold locally at each point and the renormalization yields the coefficients a , u , and 
K dependent on r and h(r) [22]. However, £ < r is violated at small r for £ > R, 
where the gradient term in (3.35) (oc V 2 ip) is indispensable. The linear relation 
tp = r£(7 3 r _7 /i(r) given in (3.34) holds for h(r) < r 136 and r > £, where the former 
condition is rewritten as r > £ + oAo^ 4 r _/3<5 / 4 . As r — > 0, the profile ip(r) becomes very 
complicated depending on A , £/i?, and the boundary condition at r = R. Detailed 
discussions will appear shortly. 
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Figure 3: Two scaling functions corresponding to form birefringence and dichroism [8]. 
The k is the wave number of probing light and £ is the correlation length. 



3.5 Electric birefringence and dichroism 

The anisotropy of the structure factor in (3.17) has not yet been measured in near- 
critical fluids, but it gives rise to critical anomaly in electric birefringence (Kerr 
effect) [16-19] and dichroism. These effects can be sensitively detected even in the 
weak field regime t c t and even for not large E\ using high-sensitivity optical 
techniques. We assume that a laser beam with optical frequency u is passing through 
a near-critical fluid along the x axis, while an electric field Eq is applied along the 
z axis. In (1.5) we have Ae yz = and the difference An = (Ae zz — Ae yy )/2n(uj) is 
written as [23, 24] 

^op 



An 



-J(k-q) 



(3.36) 



2n(u;) Jq q 2 — k 2 — iO 
where k = u/n{u) is the laser wave number in the fluid and 

A op = {de{u)/d^f/e{u). (3.37) 

Here A op is different from the static coefficient A s in (3.14). In fact A op <C A s for 
polar fluid mixtures where e(0) e{uS) at optical frequency u [40]. When r 3> r e , 
substitution of (3.18) into (3.36) gives the steady state result, 

A n 

i, 



An 



327r 2 n(u; 



F^{ki)+iFj{ki) 



(3.38) 



where the two scaling functions are plotted in Fig. 3 and are given by 
r H (r) — (I -h — ) ( A + 2 x)tan- 1 (2x)- 
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Figure 4: Normalized birefringence ReAra(i)/ReAn(0) and dichroism ImAn(t)/ImAn(0) 
vs t/tg on semilogarithmic scales for = 0.1,0.5, 1, and 2. 



For i < 1 we have -Fr(x) = Att/15 and -Fi(x) = 7rx 3 /3. In the long wavelength 
limit (A£ < 1), we find Ae cS /n(uj)An = 3(A S - 5e 2 /Q)/A op from (3.29) and (3.38), 
so AcT e fj and lim^_, An should have the same critical behavior. Experimentally, 
however, An oc r~ Zop with z op = 0.65 — 0.85 > z non was obtained [14, 18]. 

In transient electric birefringence, applied electric field is switched off (at t = 0) 
and the subsequent relaxation An(t) (for t > 0) is measured. For near-critical binary 
mixtures, this experiment has been carried out by applying a rectangular pulse of 
electric field [16-19]. If the pulse duration time At is much longer than the relaxation 
time t% of the critical fluctuations, the fluid can reach a steady state while the field 
is applied. Here t% = 6irf]^ 3 /k-BT c with fj being the shear viscosity. A remarkable 
finding is that An(t) follows a stretched exponential relaxation at short times. We 
hereafter explain our theory for weak field r c ^ t [8]. 

We assume that the relaxation of the structure factor obeys 

I{q, t) = I (q) - h{qfg e q 2 z exp(-2r,t) + • ■ • , (3.40) 

where Y q = t^ 1 _K" (g^) is the relaxation rate with K (x) = (3/4)[l + x 2 + (x 3 — 
x _1 )tan _1 x] being the Kawasaki function [22]. Thus T q = Dq 2 for g£ 1 with 
the diffusion constant D = k-QT c /Qn f)£. Here we neglect the term proportional to r c 
in (3.18) because it does not contribute to An. We substitute (3.40) into (3.36) to 
obtain An(t) as a function of k£ and t/t^, where the initial value An(0) is given by 
(3.38). In Fig. 4 we plot its normalized real and imaginary parts for various k£. The 
imaginary part decays slower than the real part. In the limit k£ — > 0, it becomes 
real and is of the form, 

An(t) 4 f°° v 2 

&A^0T G ^ ) = a " VjT ^eM-2My)t/t ( }, (3.41) 
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t/t I 

Figure 5: Comparison of the theoretical decay function (bold line) defined by (3.41) and 
data (dots) taken from Ref.[18]. The stretched exponential function exp(— 2.3a; 1 / 3 ) (dotted 
line) is a good approximation in the initial relaxation. 



where the scaling function G(x) behaves as 

G{x) 2* 1 - 2.3x 1/3 S exp(-2.3x 1/3 ) (x < 1) 

0.2x~ 3/2 (x > 1). (3.42) 

In Fig. 5 data on butoxyetbaranol + water [18] are compared with (3.41) for t < 
The agreement is excellent. It demonstrates that the theoretical G(t) is nearly 
stretched-exponential for G(t) >0.1 or for t/t^ <0.5. As another theory, Piazza et 
al. [17] derived a stretched exponential decay of An(t) on the basis of a phenomeno- 
logical picture on the distribution of large clusters. 

Let us consider the case k^ <C 1, where (3.41) is a good approximation for 
* <1/Dk 2 . At t ~ l/Dk 2 , G(t/t() becomes a very small number of order 0.2(fc£) 3 . 
In the later time region t ~^{Dk 2 ) , (3.41) cannot be used and the fluctuations with 
wave numbers of order (Dt)~ l l 2 give rise to the following birefringence signal, 

^^-m/kr 3/2 (D k H)-\ (3.43) 

Data of birefringence relaxation in Ref.[19] indicated that the decay becomes faster 
than predicted by (3.41) at long times. The same tendency can be seen in Fig. 4. In 
future such data should be compared with the theoretical curves for finite k£. We 
note that transient electric dichroism has not yet been measured for near-critical 
fluids, but was measured for a polymer solution[27]. 
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3.6 Interface instability induced by electric field 

An interface between two immiscible fluids becomes unstable against surface un- 
dulations in perpendicularly applied electric field. This is because the electrostatic 
energy is higher for perpendicular field than for parallel field. In classic papers [37] 
the instability on an interface between conducting and nonconducting fluids (such 
as Hg and air) was treated. In helium systems, where an interface can be charged 
with ions, the critical field is much decreased and intriguing surface patterns have 
been observed [38]. 

We will derive this instability in the simplest manner [39]. Let a planar interface 
be placed at z = Zq in a near-critical fluid without ions in weak field, r < and 
|t| ^> t c . If the interface is displaced by a small height ((x, y), we may set dip/dz = 
Aip5(z — z — () in (3.13), where A^ is the order parameter difference between 
the two phases. This means that the effective surface charge density is given by 
— Ei Aip E . If z is far from the boundaries, we obtain 

F dip = ^C c j dr ± J dr' ± [\C(r±) ~ CK)| 2 + |r ± - r'J 2 ]~ 1/2 , (3.44) 

where r±_ = (x,y), r' ± = (x',y'), and 

C c = A S E 2 (A^) 2 /An. (3.45) 

We expand the integrand of (3.44) in powers of [((rj_)—((rj_')] 2 . The first correction 
is negative and bilinear in (. In terms of the Fourier transformation we have 

AF dip = -±-C j dr ± J dr' ± ^ -L-^ [((r ± ) - CK)] 2 

= -^e/^|C fc | 2 , (3.46) 

where k is the two-dimensional wave vector with k — \k\ and = (2ir)~ 2 J dk. We 

have used the formula Jdr±[l — exp(ik ■ r±)]/r 3 / 2 = 4irk. Including gravity we write 
the total free energy change due to the surface deformation as 

AF surfacc = l - J k [ak 2 + gAp - C c k] | ( k | 2 , (3.47) 

where a is the surface tension, g is the gravitational acceleration, and Ap is the 
mass density difference between the two phases. The coefficient in front of | | 2 
is minimum at k — (C c /2a) l l 2 and an instability is triggered for C c > 2(agAp) l l 2 . 
For polar mixtures this criterion is typically E > \t\ u / 2 ~ 3i3 / 4 ]zV /cm on earth. 

4 Near-critical fluids with ions 

We will discuss the effects of ions doped in binary mixtures. It has long been known 
that even a small fraction of ions (salt) with c< 1 dramatically changes the liquid- 
liquid phase behavior in polar binary mixtures, where c is the mass or mole fraction 
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of ions. For small c, the UCST coexistence curve shifts upward as 

(AT) C = A lC + 0(c 2 ) (4.1) 

with large positive coefficient Ai, expanding the region of demixing. For example, 
A±/T c ~ 10 with T c ~ 300K when NaCl was added to cyclohexane+ methyl alcohol 
[41] and to triethylamine+water [42]. Similar large impurity effects were observed 
when water was added to methanol+cyclohexane [43] . In some polar mixtures, even 
if they are miscible at all T at atmosphere pressure without salt, addition of a small 
amount of salt gives rise to reentrant phase separation behavior [44, 45] . 

We consider two species of ions with charges, Ze and — e, at very low densities, 
n\(r) and n 2 (r), in a near-critical binary mixture. The average densities are written 
as (m) s = n and (712)3 — %m where (• • -) s denotes taking the space average. The 
charge neutrality condition yields Z(m) s = (n 2 ) s . The Debye wave number k and 
the Bjerrum length £ B (= 7 A for water at 300K) are defined by 

k = [Att{Z 2 + Z)ne 2 /ek B T} 1 / 2 , i B = e 2 /ek B T. (4.2) 



4.1 Ginzburg-Landau theory 

We set p = xi^E from the beginning. Then the electric potential <fi satisfies 
— V • eV(p = 4:irp with e — 1 + 4ttx and p = e(Zn 1 — n 2 ). We assume that the free 
energy F in (2.2) in the fixed charge condition is of the form, 

F = F {ip} + k B T J dr m him +n 2 lnn 2 + (wini + w 2 n 2 )ip dreE 2 , (4.3) 

where the first term depends only on ip and is given by (3.2), and the terms propor- 
tional to w\ and w 2 arise from an energy decrease due to microscopic polarization 
of the fluid around individual ions. In the neighborhood of an ion of species a 
(a = 1,2), the electric field and the polarization are given by E ind = —V(Z a e/er) 
and p ind = xEmd in terms of the local x an d e, where we write Z 1 — Z and Z 2 = — 1. 
The resultant electrostatic energy density is localized near the ion (oc r" 4 ) and its 
space integral is / dreE 2 nd /8n = Z 2 e 2 /2er a , where r a is the lower cutoff represent- 
ing an effective radius of an ion of species a. Here the screening length is assumed 
to be much longer than r a . Due to the polarization, the decrease of the electrostatic 
energy (solvation free energy per ion) is given by the Born formula[40, 46], 

AE a = (l/e-l)Z 2 a e 2 /2r a . (4.4) 

Neglecting ^-dependence of r a we estimate w a as 

w a = (A; B T)- 1 A(A J E a ) = - ei Z^ B /2e r a . (4.5) 
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The last term in (4.3) is the electrostatic free energy AF arising from slowly-varying 
fluctuations and depends on the boundary condition. As a clear illustration, if all 
the quantities are functions of z only in the charge-fixed condition, we have 

AF = 2nS [ L dz[a cx - w(z)] 2 /e(z), (4.6) 
Jo 

where a ex is the capacitor charge density and w(z) = Jq dz'p(z'). Though neglected 
here, electrostriction should also be investigated around charged particles (see Sub- 
section 3.4), which can well produce a shift of the critical composition. 

In this model the chemical potentials /i a = 5F/5n a of the ions are given by 

fi a = k B T(lnn a + 1 + w a ip) + Z a e<p. (4.7) 

In equilibrium /i a become homogeneous, leading to ion distributions, 

n a (r) = n a0 exp[-w a ij(r) - Z a e<p(r)/k B T] 

= n aQ [l-w a ^{r)-Z a e(t){r)/k B T]. (4.8) 

The second line holds when the exponents in the first line are small. Let us use the 
second line to derive Debye-Hiickel-type relations. Then n a0 = (n a ) B [l + w a {ip) s + 
Z a e{4>) s /k B T]. The potential satisfies 

e (K 2 - V 2 )0 = eiV^ • V0 + 4ire[Zn 10 - n 20 - Zn(w l - w 2 )ip}. (4.9) 

The right hand side consists of new contributions dependent on ip, where the first 
term is important in the presence of applied field. We here notice the following. 

(i) In applied electric field we can examine the ions distribution accumulated near 
the boundaries using (2.4) and (4.8) (or (4.9)). 

(ii) When w\ ^ charge distributions arise around domains or wetting layers. Let 
the interface thickness (~ £) be shorter than the Debye length kT 1 . For a spherical 
domain with radius R, for example, the electric potential 0(r) is a function of the 
distance r from the center of the sphere and then 

A(f> = (f)(0) - 0(oo) = [1 - e- KR (l + kR)]^. (4.10) 

For kR 3> 1, A0 saturates into 

$sat = -[k B T/e(Z + l)]( Wl - w 2 )A^, (4.11) 

where A*0 is the difference of ip between the two phases. For kR « 1 we find 
A0 = (/ti?) 2( l> sat . The same potential difference <3> sat arises across a planar interface. 
For example, if W\ — w 2 ~ 3, Z = 1, and Aip ~ 10 _1 , we have |$ sa t| ~ 10~ 2 V. 

(iii) In slow relaxation of Sip, the deviations 5fi a = k B T(5n a / (n a ) s + w a 5ip) + Z a e5<p 
should quickly relax to zero. Then the deviations 5n a are written in terms of Sip as 

5rii + 5n 2 = —n(wi + Zw 2 )5ip, (4-12) 



(n 2 - V 2 ) P = 



Zenivjx - w 2 )V 2 + -^-£i/t 2 E c 

4:71 



Sip, (4.13) 



where we allow the presence of electric field E to derive (4.22) below. 
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4.2 Ion effects on phase transition 

We consider small fluctuations in a one-phase state without electric field (E = 0). 
The fluctuation contributions to F in the bilinear order are written as 

SF f [1 , „ os, , ,o -sr^ / '\n n n\ 2 A 2lT 



Iq [-2 {a ° T + C ° g2)l ^' 2 + a § 2 (fe + + ^T^ 1 ^ 

1 /" [ ^ 2 + Ztt^) 2 Zn{w\ — w 2 ) 2 

2 in 



a r + C g 



1 + Z (l + Z)(g 2 + « 2 ) 



l^gl- (4-14) 



In the first line n a q and pq are the Fourier transformations of Sn a = n a — (n a ) s and 
the charge density p = e{Zn\ — n 2 ). In the second line we have expressed n a q in 
terms of ipq using (4.12) and (4.13) and minimized the first line. We introduce a 
parameter, 

7 P = ki-^2|/[(l + ^)(47rCo£ B ) 1/2 ]. (4.15) 

This number is independent of the ion density and represents the strength of asym- 
metry in the ion-induced polarization between the two components. The structure 
factor at E = in the mean field theory is written as 

1/I (q) = a (r - r ion ) + C q 2 [l - ^ / '(q 2 + k% (4.16) 

where 

Tion = K + Zw 2 fn/{1 + Z)a . (4.17) 

We draw the following conclusions. 

(i) If 7 P < 1, Io(q) is maximum at q = and the critical temperature shift due to 
ions is given by T c T ion in the form of (4.1). As a rough estimate for Z — 1, we set 
0-0 X n ~ £+ ^ ~ c i where c is the mass or mole fraction. Then r ion ~ (wi + w 2 ) 2 c. 
If \wi + w 2 \ ~ 3, this result is consistent with the experiments [41, 42]. In future 
experiments let l/C I (q) vs q 2 be plotted; then, the slope is 1 — j 2 for q k and 
is 1 for q ^> k. This changeover is detectable unless 7 P 1. 

(ii) The case 7 P = 1 corresponds to a so-called Lifshitz point [47], where 1/I(q) — 
1/1(0) ocq 4 /(q 2 + K 2 ). 

(iii) If 7 P > 1, the structure factor attains a maximum at an intermediate wave 
number q m given by 

Qm = (7 P - 1) 1/2 ^ oc n 1/2 . (4.18) 
The maximum structure factor Io(q m ) diverges as r — > r/ on , where 

< n = ^ion + (7 P -l) 2 2 , (4-19) 

with the aid of C = a £+ . A charge- density- wave phase should be realized for 
t < r/ on . It is remarkable that this mesoscopic phase appears however small n is 
(as long as 7 P > 1 and q m L 3> 1, L being the system length). Here relevant is the 
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coupling of the order parameter and the charge density in the form oc ipp in the 
free energy density, which generally exists in ionic systems. This possibility of a 
mesoscopic phase was already predicted for electrolytes [48], but has not yet been 
confirmed in experiments. 

Furthermore, we note that 7 P increases with increasing Z. As an extreme case, 
we may add charged colloid particles with Z ^> 1 and radius R in a near-critical 
polar mixture, where W\ grows as Z 2 from (4.5) and 

7p = (£B/lQ7rC ) 1 / 2 (e 1 /e )Z/R. (4.20) 

Note that Z£ +0 /R can be made considerably larger than 1 with increasing R. Notice 
that the ionizable points on the surface is proportional to the surface area 4nR 2 . 

4.3 Nonlinear effects under electric field 

In most of the previous experiments, a pulse of strong electric field has been applied. 
For example, the field strength was E = 10 4 V/cm and the pulse duration time was 
At = 1ms at T ~ 300K [13, 15]. Let us set Z — 1. As can be seen from the general 
relation (2.18), we can neglect ion accumulation at the capacitor plates if 

a cx = eEq/Att > env dri At or At < t dri . (4.21) 

Here v dri = £ _1 e-E is the drift velocity, ( is the friction coefficient related to the 
diffusion constant by = k^T/(, and t dT \ — {k 2 D{)~ 1 is the drift time. For At > 
tdri the electric field far from the boundaries are shielded if a cx < enL. Typical 
experimental values are A ~ 10~ 6 cm 2 /s, ( ~ lCT 8 g/s, e ~ 50, and E = 10 4 V/cm, 
leading to v<iri ~ 1 cm/s. Then the condition (4.21) becomes nAt <^ 10 12 cm~ 3 s and 
can well hold in experiments [49]. Under this condition and in a one-phase state, 
the bulk region remains homogeneous and the electric field is not yet shielded. If we 
are interested in slow motions of ip, we may assume (4.12) and (4.13) to obtain the 
decay rate of ipq in the form Tq = L q 2 /I(q) (L being the kinetic coefficient) with 

HQ)- 1 = Hq)- 1 + [S2 + A s q 2 z /(q 2 + k 2 )] (E 2 /4rrk B T), (4.22) 

where Io(q) is given by (4.16). If the pulse duration time At is much longer than the 
relaxation time 1/Tg, the structure factor is given by I(q). Here q 2 = q 2 z /q 2 in (3.17) 
is replaced by q 2 /(q 2 + k 2 ) in the presence of ions due to the Debye screening. Let 
us consider the nonlinear dielectric constant in (3.29) and the electric birefringence 
lim^^o An in (3.36) on the order of E 2 . For < 1 the ion effect is small, but for 
k(, > 1 they should behave as 

Ae cS /A s ~ An/A op oc E 2 k~\ (4.23) 

This crossover occurs for c ~ n£ij_ > r 2v in the weak field regime r 3> t c . 
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We also comment on the Joule heating. While the ions are drifting, the temper- 
ature increasing rate is given by 

CjT = 2nv Ali eE = 2ne 2 E 2 /( = eE 2 /4irt drh (4.24) 

where C is the specific heat per unit volume at constant volume and composition. 
For At > t d ri the temperature will increase by (AT)j oule ~ eEq/AttC. By setting 
C ~ fcutn), we obtain dr/dt ~ 10" 19 n ~ lO" 6 ^ s" 1 at E = 10 4 V/cm. The heating 
is negligible for very small nAt or for not very small r. 

As other nonlinear problems involving ions, we mention transient relaxation of 
the charge distribution after application of dc field, response to oscillating field, and 
effects of charges on wetting layers and interfaces between the two phases. 



5 Liquid crystals in electric field 

In liquid crystals near the isotropic-nematic transition, the order parameter is the 
symmetric, traceless, orientation tensor ip = {Qij} (which should be distinguished 
from the electric charge Q on a capacitor plate). The polarization and the electric in- 
duction are written as p —X E and D =e -E, respectively, where the polarizability 
tensor X depends on and the local dielectric tensor e = 1 + 4ir X reads 

where Q 2 = Q% hereafter. In the nematic phase we may set Qij = S(riinj — 5ij/3) 
in terms of the amplitude S and the director n. Then, 

Eij = e' (S)Sij + e'^SfyiUj, (5.2) 

where e' = Eq + S ,2 (2e 2 i/3 + £22/9) and e[ = e±S + E22S 2 /3. These relations are 
analogous to (1.2) for near-critical fluids. An important difference is that the tensor 
{Qij} is n °t conserved and its average is sensitive to applied electric field, while the 
average order parameter in near-critical fluids is fixed or conserved. 



5.1 Pretransitional growth 

Here we examine the effect of the field-induced dipolar interaction. If we assume 
no mobile charges inside the fluid, the free energy functional is given by F = E + 
J drE ■ D/8-K at fixed capacitor charge, where 

F =f dri^p-Q 2 + Qa + y YpkQiif + y E ViQijVkQkj + • • •] (5.3) 

j ijk ijk 
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is the Ginzburg-Landau free energy for {Qij} and A is a space-dependent Legendre 
multiplier ensuring J2jQjj — [20]. For weakly first order phase transition, the 
coefficient A(T) — a\(T — To) becomes small as T approaches To. In To the higher 
order terms are not written explicitly. Analogously to (3.9) we find 

5F = 5F - J drE ■ 5s ■ E/8n. (5.4) 

Thus, 

-F=[A(T)Q ij + \5 ij + ---]--±E i E j , (5.5) 



SQij L v lJ 1 8tt 

where the gradient term is omitted. At fixed capacitor potential, on the other hand, 
the appropriate free energy is G = T — / drE ■ D/8n but 5G/5Qij is of the same 
form as 5F/5Qij in (5.5). In equilibrium disordered states with T > T , we set 
{SF/SQij) = to obtain 

(Qij) = ^(E 0i E 0j - -E 2 ^)-^ + • • • , (5-6) 

where E is the average electric field assumed to be along the z axis. On the other 
hand, analogously to (3.27), the macroscopic static dielectric constant is given by 

e e a = (e zz ) - f -^(| Y>9j £ M\ 2 ), (5-7) 

where Sij{q) are the Fourier transformation of £ij(r). As T — > T , the average (e zz ) 
yields the dominant contribution to the nonlinear dielectric constant [20], 

Ae cff = e\El/l2nA{T), (5.8) 

in agreement with the experiments [21]. The contribution from the second term in 
(5.7) is smaller than that in (5.8) by (e 1 / 'e ) 2 k B T 'A(T) 1 / 2 / 'L^ 2 for L 2 ~ L x . 
The electrostatic free energy AT = T — T up to of order Q 2 is written as 

AF = ~§J dre zz (Q)+A s ^ /^|£^(<?)| 2 , (5.9) 

where A s = e 2 /e and the dipolar interaction, the second term, is expressed in 
terms of the Fourier transformations Qij(q). The correlation functions of Qij(q) 
in disordered states depend on the direction q~ x q even in the limit q — > 0. For 
simplicity, for q z = we obtain 



k T 

\im{Q iz (q)Q jz (q)*) = — 



q l r + r c 



(hj=x,y), (5.10) 



where r = 2A(T) - (2e 2 i + e 2 2)E$/47t and r c = A s E%/4n. 
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5.2 Director fluctuations in nematic states 



We consider nematic states considerably below the transition, where we may neglect 
the fluctuations of the amplitude S [20]. Hereafter we rewrite e' and e[ in (5.2) as 
Eq and E\ for simplicity. If e is positive, the average orientation of the director can 
be along the z axis from minimization of the first term in (5.8). Then the deviation 
Sn perpendicular to the z axis undergo large fluctuations at small wave numbers. 
In electric field AF in (5.9) becomes 



AF = ^El f L\Sn(q)\ 2 + A' s \\q ■ 5n(q)f 

07T J Q L Q 



1 - ^ £l \5n( q )\> + A' s ± 
where A' s = A S S 2 . For general q we obtain the correlation functions, 



(5.11) 



GiM) = (Sn^Snjiq)*) = k B T^ + k B T^(- - -) (ij = x,y), (5.12) 

r 2 qj_ Vri r 2 J 

where q\ = q\ + q 2 , and 

r^Ksqj + ^ql + j-^ + A'^El r 2 = K 3 q 2 z + K 2 q\ + j-e.E 2 . (5.13) 

The Ki, K 2 , and K 3 are the Frank constants. If E\ >A' S , the correlation length £ is 
of the following order, 

f ~ (AtcK/s^E^ 1 , (5.14) 

where K represents the magnitude of the Frank constants and k B T / K is a micro- 
scopic length. The scattered light intensity is proportional to the following [20] , 



\f-7 (q)-i\ 2 ) = k B Te 2 1 



|a|7r 2 + (l/r 1 -l/r 2 )|q-a| 2 /g 2 



(5.15) 



where * and / represent the initial and final polarizations. The vector a is defined 
by 

&x = izfx + fzix-, Q>y = izfy + fz^yi a z = 0. (5.16) 

If A' s > 0, the intensity depends on q^q even in the limit q — > 0. 

In nematic states the average (e^) is anisotropic, leading to large intrinsic bire- 
fringence. We here consider form dichroism arising from the director fluctuations, 
where the laser wave vector k is along the x axis and the average director is along the 
z axis. We assume the relation = e 5ij + e-jiirij at optical frequencies using the 
same notation as in (5.2). As a generalization of (1.5), the fluctuation contribution 
to the dielectric tensor for the electromagnetic waves is of the form, 



e o Jq q k ~ *0 



(5.17) 
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where [• • -]± denotes taking the tensor part perpendicular to k. This expression 
reduces to (1.5) if e is diagonal and (1.2) holds. Here Ae yz = and the imaginary 
part of An/n = (Ae zz — Ae yy )/2e becomes 



ImAn k 3 



n 



32tt 2 



e l /e ) 2 Jd^G xx (£) + q 2 z G yy (£)- £ q^GM 



i,j=x,y 



(5.18) 



where £ = k — kq and / dQ ■ ■ ■ denotes integration over the direction q = q l q. 
Using (5.12) we can make the following order estimations, 

— ~ io- 2 ( £l / £o ) 2 (fer/ir)^ 3 (^«i) 

n 

~ 10-\e 1 /e ) 2 (k B T/K)k (A£ >1). (5.19) 

The form dichroism here is much larger than that in (3.38) for near-critical fluids. 
As far as I am aware, there was one attempt to measure anisotropy of the turbidity 
in an oriented nematic state [50]. 

5.3 Orientation around a charged particle 

As another example, we place a charged particle with radius R and charge Ze in a 
nematic state, where n is aligned along the z axis or n — > e z far from the particle. 
Let the density of such charged particles be very low and its Coulomb potential 
be not screened over a long distance A (which is the Debye screening length if low 
concentration salt is doped). From (5.4) the free energy change due to the orientation 
change 5n is given by 



j (h 



KV 2 n + ^-(E-n)E 

47T 



Sn, (5.20) 



where we assume the single Frank constant K (K 1 = K 2 = K 3 = K). If the 
coefficient E\ is considerably smaller than £ 0) the electric field E near the particle 
is of the form — V{Ze/ e^r). Then, for E\ > (or E\ < 0), n tends to be parallel 
(or perpendicular) to r — r~ 1 r near the charged particle. We assume that n is 
appreciably distorted from e z in the space region R < r < A. For e 1 > the 
decrease of the electrostatic energy is estimated as 



2e 2 \R 

analogously to (4.4). For e\ < 0, e\ in (5.21) should be replaced by |ei|/3 (because 
the angle average of (e z ■ r) 2 is 1/3). The Frank free energy is estimated as 

F ~ irKe. (5.22) 
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Figure 6: The director field in two dimensions for £i/eq = 0.4 in (a) and for E\/eq = —0.4 
in (b) around a charged wire at the origin. It tends to be along the y axis (in the vertical 
direction) far from the origin. The spacing between the adjacent bars is 2^b- 



We determine I by minimizing Fq + AF to obtain 

£=(\e 1 \/2TTs 2 K) 1 / 2 Ze. 



(5.23) 



The condition of strong orientation deformation is given by R < I or 

R/Z < (Is^neo^ie^T/K) 1 / 2 . (5.24) 

If this does not hold, the distortion of n becomes weak. 

Fig. 6 illustrates the deformation of n in two dimensions. We have numerically 
solved 



n x 



KV 2 n 



— (E ■ n)E 

4vr ; 



0. 



(5.25) 



under —V- e -V0 = Airp by assuming n = (cos 9, sin 0) (or n z = 0). A charge is 
placed in the hard-core region (x 2 + y 2 ) 1 ^ 2 < R. In three dimensions this is the 
case of an infinitely long charged wire with radius R and charge density cr, in which 
all the quantities depend only on x and y. The solution can be characterized by 
the three normalized quantities, £i/£q, cr* = a / \sqK) 1 I 2 , and R/£b- Here we set 
a* = 2.4 and R/£b — 2. We discretize the space into a 200 x 200 lattice in units of 
£b under the periodic boundary condition in the x direction, so the system width is 
L = 200^b- The electric potential vanishes at y = and L. 



6 Concluding remarks 

(i) The Ginzburg-Landau theory in Section 2 generally describes how the electro- 
static interactions arise depending on the boundary condition (in the presence or 
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absence of the capacitor plates). It can be used to investigate electric field effects 
at various kinds of phase transitions in fluids and solids. 

(ii) A brief review has been given on the dielectric properties and the electric field 
effects in near-critical fluids and liquid crystals. The Debye-Kleboth experiment on 
the critical temperature shift was performed many years ago and they neglected 
the dipolar interaction in their theoretical interpretation. As regards the nonlinear 
dielectric constant Ae c s and the birefringence An, we cannot explain their experi- 
mental exponents, z non ~ 0.4 and z op = 0.65 — 0.8, whereas the common exponent 
v{\ — 2i]) = 0.58 has been predicted for them. To resolve these issues, scattering 
experiments to check the anisotropic structure factor (3.f7) are most needed. 

(iii) New predictions have also been made on the critical temperature shift due to 
electric field, the nonlinear dielectric constant, the ion effects in binary mixtures, 
and the fluctuation intensities and the form dichroism in liquid crystals. In particu- 
lar, in near-critical polar mixtures with ions, we have examined charge distributions 
and potential differences around two-phase interfaces, the critical temperature shift 
due to ions, and the scattering intensity. The condition for a charge-density-wave 
phase has been examined for general multivalent ions. 

(iv) Effects of oscillating electric field are also worth studying particularly for ionic 
systems. Appreciable critical anomaly can be seen in the frequency-dependence of 
the dielectric constant [49]. Large dynamic electric birefringence was observed in 
polyelectrolyte solutions [51]. By this method we can neglect accumulation of ions 
at the boundaries (but the Joule heating may not be negligible). 

(v) In Section 5 we have examined the deformation of the nematic order around a 
charged particle. The charge- induced orientation is intensified with decreasing the 
particle radius R and/or increasing the charge number Z. This is in marked con- 
trast to the surface anchoring of a neutral particle [52, 53], which can be achieved 
for large radius because the penalty of the Frank free energy needs to be small. For 
charged colloidal suspensions, furthermore, the counterions themselves can induce 
large deformation of the nematic order (because of their small size) and tend to 
accumulate near the large particles. These aspects should be examined in future. 

(vi) Stronger electric field effects have been observed in polymeric systems than in 
low-molecular-weight fluids. As such effects, we mention field-induced anisotropy in 
light scattering from a polymer solution [32], lamellar alignment in diblock copoly- 
mers [54], and large dielectric response in a surfactant sponge phase [55]. 
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